Conductance properties of nanotubes coupled to superconducting leads: signatures of 

Andreev states dynamics 
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We present a combined experimental and theoretical analysis of the low bias conductance prop- 
erties of carbon nanotubes coupled to superconducting leads. In the Kondo regime the conductance 
exhibits a zero bias peak which can be several times larger than the unitary limit in the normal 
case. This zero bias peak can be understood by analyzing the dynamics of the subgap Andreev 
states under an applied bias voltage. It is shown that the existence of a linear regime is linked to 
the presence of a finite relaxation rate within the system. The theory provides a good fitting of the 
experimental results. 

PACS numbers: PACS numbers: 73.63.-b, 73. 21. La, 74.50.+r 



Carbon nanotubes connected to metallic electrodes al- 
low to study many different regimes in mesoscopic elec- 
tron transport exhibiting either ballistic or diffusive be- 
havior [l|,|2 [; Luttinger liquid features J3[, etc. Several ex- 
periments J| have demonstrated that carbon nanotubes 
weakly coupled to normal leads can also exhibit Coulomb 
blockade and Kondo effect, i.e. the characteristic physics 
of quantum dots (QD). Very recently the research on the 
transport properties of this system has been extended 
with the inclusion of superconductivity on the leads p| . 
In particular Buitelaar et al. ^ have achieved a physical 
realization of a S-QD-S system using carbon nanotubes 
coupled to Al/Au leads. Their results, showing an en- 
hancement of the conductance in the Kondo regime with 
respect to the normal case, have attracted considerable 
theoretical attention 0, H ■ 

From a theoretical point of view transport in a volt- 
age biased S-QD-S system poses a challenging problem 
due to the interplay between electron correlation effects 
and multiple Andreev reflection (MAR) processes which 
provide the main mechanism for quasiparticle transport 
at small bias. This problem has been addressed by sev- 
eral authors at different levels of approximation. Thus, 
in Refs. [llj numerical results were presented for the S- 
QD-S system neglecting electron correlation effects. This 
approach has also been used by Buitelaar et al. [7| for 
analyzing the experimental results for the subgap struc- 
ture of Ref. la]. On the other hand, Avishai et al. |lO| 
have introduced the effect of electron interactions at the 
level of the slave boson mean field approach and obtained 
numerical results for the current and the noise spectrum 
in the Kondo regime. More recently, some of us have 
shown that the conductance properties of a S-QD-S sys- 
tem in the Kondo regime can be understood in terms of 
the dynamics of subgap states (Andreev states) under an 
applied bias voltage Q. 

In the present work we combine experimental and theo- 
retical analysis in order to identify signatures of Andreev 



FIG. 1: Typical 2-terminal device geometry. For the mea- 
surements presented here, the electrode spacing is 0.25 fim 
and the MWNT length 1.5 ^m. The (oxidized) Si substrate 
is used as a gate electrode. 



states dynamics in the low bias conductance of carbon 
nanotubes coupled to superconducting leads. We shall 
first describe the setup used for the conductance mea- 
surements and discuss the main experimental results. We 
then introduce the theoretical model in which we assume 
that the system can be appropriately analyzed as a single 
level quantum dot. In contrast to Ref. |^, in the present 
calculations we explicitly take into account the presence 
of a finite inelastic relaxation rate within the system. As 
will be shown the low bias conductance properties are ex- 
tremely sensitive to this rate. The comparison with the 
experimental results allows to determine the size of this 
parameter in the actual system as well as to confirm the 
main predictions of the theory. 

The device we consider consists of an individual mul- 
tiwall carbon nanotube (MWNT) of 1.5 /im length be- 
tween source and drain electrodes that are separated by 
250 nm, Fig^ The lithographically defined leads were 
evaporated over the MWNT, 45 nm of Au followed by 
135 nm of Al. The degenerately doped Si substrate was 
used as a gate electrode. 

The device is first characterized with the contacts 
driven normal by applying a small magnetic field of 
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FIG. 2: Gray scale representation of the differential conduc- 
tance dl /dV versus source-drain (Vsd) and gate voltage (Vg) 
with the leads in the normal state (darker — more conduc- 
tive). The dashed white lines outline the Coulomb blockade 
diamonds and the Kondo ridge. The black dashed line trace 
is the dl /dV in the middle of the Kondo ridge (position B) , 
from which Tk can be estimated. The line trace in the bottom 
graph is the linear-response conductance. 



B — 26 niT. The critical field of the electrodes was ex- 
perimentally determined to be ~ 12 niT. Figure |21 shows 
a gray scale representation of the measured differential 
conductance dl/dVsd versus source-drain [Vsd) and gate 
voltage (V^) at T = 50 mK with the leads in the normal 
state. Measured over a large range of Vg we observe a 
series of truncated low-conduction Coulomb blockade di- 
amonds linked by narrow ridges of high conductance of 
which Fig[5]is a small part, see also Ref._6j. From these 
and other Coulomb diamonds at different Vg we obtain 
an average charging energy Uc = 0.45 meV and an aver- 
age level spacing 5e = 0.45 meV, with Se varying between 
0.3 and 0.7 meV. 

The high conductance ridge around Vsd = mV in 
Fig|21 is a manifestation of the Kondo effect occurring 
when the number of electrons on the dot is odd and the 
total spin S = 1/2. At T = 50 mK the Kondo state 
seems to be fully developed and has reached the unitary 
limit with a conductance maximum of 2e^ jh. The Kondo 
temperature Tk can be estimated from the width of the 
Kondo ridge out of equilibrium [l2|. The full width at 
half maximum (FWHM) corresponds to ^ fesTx which 
yields a Kondo temperature for the ridge in Fig[21of about 
1.86 K. In total we have observed 12 Kondo ridges with 
Kondo temperatures varying between 0.52 and 3.34 K. 

When the magnetic field is off, the leads are super- 
conducting and the conductance pattern is dramatically 
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FIG. 3: Differential conductance dl /dV in the middle of a 
Coulomb diamond (left) and in the middle of the Kondo ridge 
(right). These correspond to positions A and B in Fig|5| re- 
spectively. The difference between the superconducting state 
of the leads (solid line) and normal state (dashed line) is 
clearly observed. The peaks in the dl /dV in the supercon- 
ducting state are due to MAR. 



different. This is illustrated by Fig|31 which compares 
dl/dV in the normal (dashed lines) and superconduct- 
ing state (solid lines) at the Vg positions A and B of 
Fig|21 Whereas the dl/dV in the normal state has no 
pronounced features (except for the Kondo resonance at 
position B), a series of peaks appear in the superconduct- 
ing state. These are the signature of MAR. The peaks 
at ±0.2 mV correspond to the superconducting gap of 
2A and mark the onset of direct quasiparticle tunnel- 
ing. The peaks at lower Vsd involve one or more Andreev 
reflections. The observed subharmonic gap structure is 
discussed in more detail in Ref. 7] . The superconducting 
gap energy of A = 0.1 meV is slightly less than the ex- 
pected 0.18 meV for bulk aluminum which we attribute 
to the intermediate Au layer (necessary to obtain good 
electrical contact to the MWNT). 

Whereas the conductance at Vsd — mV in the su- 
perconducting state is less than its normal state value 
in the middle of an 5* = Coulomb diamond, the sit- 
uation is very different for the Kondo region. Here the 
conductance increases significantly and a resonance re- 
mains around Vsd = 0, albeit narrower than in the normal 
state. This has been observed for all Kondo resonances, 
provided Tk ^ A, see Ref.lQj]. As stated in the intro- 
duction, here we focus on the low-bias region of those 
Kondo resonances in the superconducting state and try 
to identify signatures of Andreev states dynamics. 

To describe theoretically a S-nanotube-S system we 
shall first assume that the level spacing Se in the nan- 
otube is large compared to the superconducting gap in 
the leads A. In the actual case of the present experi- 
ments A ~ 0.1 meV and Se ~ 0.45 meV, which means 
that this approximation is roughly correct. The problem 
can thus be mapped into a single-level Anderson model 
with a model Hamiltonian H = Hl + Hji + Ht + Hq: 
where Hlh describe the left and right leads as BCS su- 
perconductors, Ht = X]fe""-r D^ucl^do- + h.c, is the 
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term coupling the dot to the leads, d^ and cj, ^ being cre- 
ation operators for electrons in the dot and in the leads 
respectively. Hq = X^o- ^o"-o- + Un-jfi^ is the uncoupled 
dot Hamiltonian characterized by the dot level position 
eo and the Coulomb interaction U. It is assumed that 
the coupling of the dot to the leads in the normal state 
can be described by energy independent tunneling rates 

A crucial point in the analysis of this model is how to 
deal with the electron correlation effects. In Ref. |1^ a 
perturbative approach was used to study the zero- voltage 
case. It was found that the AS's in the Kondo regime 
(Tr: > A) have essentially the same phase dependence as 
in the non-interacting case but with a reduced amplitude. 
This is consistent with a Fermi liquid description of the 
normal state with renormalized parameters Eq and F^ fj, 
where F? -^Ft, = 2F* ~ Tk fixes the width of the Kondo 



resonance. Although different approximations differ in 
the way to relate these quantities to the bare parameters, 
the description in terms of the renormalized ones can be 
considered as universal [lj|. In the opposite limit {Tk < 
A), the system is well described within the Hartree-Fock 
approximation pXl Ha. Il5l| . 

This fully phase-coherent picture is limited in the ac- 
tual system by the presence of inelastic scattering mech- 
anisms not considered in the present model like electron- 
phonon interactions. We shall later discuss its effect in 
the quasiparticle current by introducing a phenomcno- 
logical inelastic relaxation rate rj. 

In the zero-bias limit the Andreev states are deter- 
mined by the poles of the dot retarded Green function 
which, in the approximation mentioned above can be 
written as 
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-(tj-t-iO+)/v/A2-(cj-f iO+)2 
- i0+)2 corresponding to the 



A/ VA2 - (u; ^ 
dimensionless Green functions of the uncoupled elec- 
trodes. In the electron-hole symmetric case (ej — and 
F^ = F|j.) the AS's are then determined by the equation 
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(2) 



The solutions of Eq. |(5J are plotted in Fig. 0J As can 
be observed the AS detach from the continuous spectrum 
and behave approximately as ujs ~ ± A cos 0/2 with A < 
A. The energy interval between the AS and the edges of 
the gap is given by A - A ~ 2A^/{2T*)^ for F* > A. 
On the other hand at the crossing between the two AS's 
(0 = ±7r) A can be approximated by A/(l -I- A/2F*). 

In a non-symmetric situation (cq ^ and/or F^ 7^ F|j) 
the two ballistic states are coupled and there appear an 
upper and a lower bound state with an internal gap of the 
order of 2Ay/T^^, where r = 4F*F^/ [{e^)^ + (F*)^] 
is the normal transmission at the Fermi energy. This is 
similar to the behavior of the AS's in a point contact with 
finite reflection probability, except for the detachment 
from the continuous spectrum (see Fig. ^J. 

In the voltage biased case it can be assumed that the 
description in terms of the renormalized parameters e|5 
and F^ ^ remains valid in the regime eV <^ A ~ T^ in 
which we are interested. Within the present model the 
current operator between the leads and the dot is given 
by 
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FIG. 4; Andreev states for a S-QD-S in the Kondo regime. 
Dashed and full lines correspond to a symmetric and non- 
symmetric situation respectively. 






^t' ^ J 1^ i^t^^'lcdr^ 



h.c 



(3) 



For calculating the average current we use the Keldysh 
formalism following previous works |lfi| . The main quan- 
tities to determined are the Keldysh Green functions in 



Nambu space given by Gf, {t,t') — i < ■ip*-{t')ipi{t) >, 

where ipi are the field operators in Nambu space defined 
as 



V'» = f ^t^ j , V'l = ( c\^ C,i ) 



The average current can then be written as 
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where t^ = t^itTz. One can integrate out the leads to 
find an expression of the current only in terms of the 
dot Keldysh Green function G J~ , which is related to the 
retarded and advanced Green functions by 
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where gj^ ^ are the Keldysh Green functions for the un- 
coupled electrodes, whose Fourier transform is given by 



L,Rt 
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{Lij)[g°'{uj) — g^{uj)\, Up [oj) being the correspond- 
ing Fermi factor. In Eq. lO integration over internal 
time arguments is implicitly assumed. 

Using the different relations between Green functions 
in Keldsyh space the current between the dot and the 
left lead can be written as 
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FIG. 5: Low bias current-voltage characteristic for the sym- 
metric case. Panel (a) corresponds to the fully phase-coherent 
case (77 = 0) and different values of r*/A = 10, 5, 4, 3, 2, 1 
from top to bottom. Panel (b) corresponds to r*/A = 2 and 
different values of r;/A as indicated. The dotted lines show 
the analytic description of the eV < rj regime given by Eq. 
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where f^ — t^GDtt- In order to evaluate the current 
under fixed bias voltage a double Fourier transformation 
of the temporal arguments in the dot Green functions is 
performed, which leads to a set o f algebr aic equations for 
the different Fourier components [ulllg- This procedure 
allows to evaluate the current as I{t, V) = [II — Ir)/"^ = 
Y.nIri{V)exp[in(j){t)], where (j)(t) = 2eVt/h + ^q. In the 
present case we shall be interested in the dc component 
h{V). 

The behavior of the dc current at low bias is shown 
in Fig. |S1 When F* » A the IV characteristics of a 
ballistic quantum point contact is gradually recovered. 
As F*/A is reduced the low bias conductance is gradually 
suppressed. When F* ~ A, MAR oscillations with a 
period A/n, where n is an integer, start to be observable 
in the IV characteristic Q- 

These results correspond to the fully phase-coherent 
case (77 = 0) and, as can be observed, exhibit a highly 
non-linear behavior where a linear regime is strictly never 
reached. The presence of a finite 77 in the actual system 
limits the observation of this behavior, giving rise to a 
linear regime in the limit V -^ Q. This is illustrated in 
the right panel of Fig. where the IV characteristic for 
F* = 2A is represented for different values of 77. As can 
be observed, the current at low bias is extremely sensitive 
to the actual value of this parameter. Notice also that the 
oscillations with period A/ti are suppressed for eV < rj. 



These results can be understood more deeply by ana- 
lyzing the dynamics of the Andreev states in the presence 
of an applied bias taking into account the finite relax- 
ation rate. Let us first consider the symmetric case. In 
the limit eV < 77 we can obtain the non-equilibrium oc- 
cupation of the Andreev states, p , using the Boltzman 
equation in the relaxation time approximation, i.e. 
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where p^^''^ = np [±^5(0)] and r = h/2r]. In the 
stationary state these quantities can be expanded as 
p ~ ^jj^ pn exp inuj^t, where loq — eV/fi and we obtain 
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At zero temperature the Fourier components of the 
equilibrium distribution are simply given by P2fc+i ~ 
(— l)'^/(7r(2fc -I- 1)). To calculate the current we first no- 
tice that within this approximation 






(10) 



Averaging this expression to obtain the dc current, as- 
suming that ws — A cos </)/2, yields 



2e2A 
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(11) 



This expression indicates that the conductance reaches 
a maximum at F = where it takes the value G'(O) = 
2e^A/(/i77) and decreases with V defining a zero-bias con- 
ductance peak with a typical width fixed by ^ 77. In panel 
(b) of Fig. [3 we compare the results of this approxima- 
tion with the full numerical ones. As can be observed 
Eq. pi(l describes correctly the behavior for eV < rj. 

This description can be straightforwardly generalized 
to a non-symmetric situation. In this case the effect of 
the transition between the lower and the upper state can 
be analyzed as a Landau-Zener process with probability 
p = exp [—TTrA/eV), where r = 1 — r |l7l |. The current 
in Eq. 1)11(1 will be then suppressed by a factor p. 

On the other hand for finite temperature T, we obtain 
/o(T) = Io{T = O)F(r), where Io{T = 0) is given by Eq. 
(fTT|l and 



F{T) = ^ V 



yr 



(12) 



where P = l/ksT and x^ = (2m + l)7r/(/3A). The ab- 
sence of any dependence on 77 or F in F{T) indicates that 
the main effect of temperature is to introduce a global 
scale factor reducing the magnitude of the conductance. 

This simple theory can be used to analyze the low bias 
conductance peak in the experimental results. As shown 
in Fig. El this peak can be well fitted by the expression 
(|ll|l using rj and A as free parameters. The values of 
77 thus obtained are rather constant (between ^ 10 and 
^ IbfieV) while A strongly depends on Tk- This be- 
havior is consistent with our theoretical model in which 
the inelastic relaxation rate is a fixed parameter set by 
some microscopic mechanism not explicitly included in 
the model. Notice that rj is considerable smaller than 
both A and Se. On the other hand the theory predicts 
an increase of A as a function of Tx roughly given by 
A/ A ~ 1/(1 + A/Tr-). In the lower right pannel of Fig. 
Elthis prediction is compared with the fitted experimental 
values. Notice that, in spite of the rather large dispersion 
in the data, both the order of magnitude and the general 
trend are in agreement with the theoretical predictions. 

In conclusion we have presented a combined theoretical 
and experimental analysis of electron transport in carbon 
nanotubes coupled to superconducting leads. The differ- 
ential conductance in the Kondo regime exhibits a zero 
bias peak which can be related to the dynamics of the 
subgap Andreev states under an applied bias voltage. A 
main ingredient that has to be included in the model in 
order to account for the experimental results is a finite re- 
laxation rate 77 which gives rise to a linear regime in the 
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FIG. 6: Differential conductance in the middle of 5 different 
Kondo ridges with decreasing Tk- The red lines are fits to 
the data using Ea llll The lower right pannel shows the fit- 
ted values of A vs Tk compared to the theoretical prediction 

A/A ~ 1/(1 -fA/TK). 



low bias limit eV < rj. The different mechanisms that 
can contribute to this finite rate can be either intrinsic 
(like electron-phonon coupling within the nanotube) or 
extrinsic (for instance due to phase diffusion like in a 
RS J model 9] ) . Further experiments could be of interest 
in order to elucidate this issue. 
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